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Abstract 

We find a parametric resonance in tlie GHz range of the DNA dynamics, 
generated by pumping liypersound . Tliere are localized phonon modes caused 
by the random structure of elastic moduli! due to the sequence of base pairs. 



1 Introduction: phonon modes of the DNA 

Molecule of the DNA has unusual elastic properties owing to its helicoidal symmetry 
and the sequence of base-pairs. From the physical view-point, i.e. neglecting its 
genetic information, the latter looks random. The vibrational dynamics of the DNA 
has an important bearing upon biological phenomenae in cell, [I], see also[5]. A 
specific feature of the dynamics is localized motions in the duplex which spread only 
over several base-pairs. Their existence has been confirmed by the experimental 
research using the Raman scattering, [3], [4], [5], [7] , the far- infrared absorption, 8^,\9^, 
and the Brillouin scattering [10]. In paper [11] the vibrational modes have been 
studied using the submillimeter-wave absorption spectroscopy in the range of 0.01 — 
10 THz. Thus, Woolard et al, [llj, have found multiple dielectric resonances in the 
long- wavelength portion of the submillimeter-wave regime, i.e. ^ 1 — 30cm~^, which 
they ascribe to phonon modes of the DNA. These results provide valuable methods 
of bio-diagnostics, [12) . 
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Theoretical study of phonon modes of the DNA had revealed elastic excitations 
of the duplex which may correspond to the approximate helicoidal symmetry of a 
molecule of the DNA, and its random elastic structure, [Hj, [TS]. It is important 
that phonon modes of the DNA are believed to be strongly attenuated owing to 
an interaction with ambient medium. The effect can be mitigated, to some extent, 
by preparing samples of appropriate character. Thus, films formed by molecules of 
the DNA appears to be less prone to the attenuation as regards phonon modes. In 
contrast, it is especially strong in liquid solutions of the DNA. It should be noted 
that the estimates are based on the Navier-Stokes hydrodynamics in sm6-GHz range, 
even more so, the classical Stokes formula for the viscous drag at small Reynolds 
number. But the above arguments are wrong in the GHz range. In fact. Van Zandt, 
[16], showed that within the framework of the Maxwell hydrodynamics, there could 
exist underdamped phonon modes. Thus, whether the phonon modes are damped or 
not so damped, remains to be seen. 

In this paper we follow the analysis performed by Chia C.Shih and S.Georghiou 
who have put forward powerful arguments in favour of the existence of underdamped 
vibrational modes of a molecule of the DNA, besides the overdamped ones P^. The 
authors visualize a molecule of the DNA as a duplex of two strands formed by sugar- 
phosphate backbones framed by base-pairs located inside the strands, and assume that 
the bases are shielded by sugar — phosphates from the bombardment by molecules of 
solvent, the base-pairs being influenced by the ambient medium indirectly, through 
their interaction with the backbone. The motion of the bases has the shape of li- 
brations inside the cages formed by the sugars of the backbone. Consequently, the 
dynamics of the backbone is damped owing to the strong attenuation caused by the 
medium, whereas that of the base-pairs turns out to be underdamped. Their angular 
frequencies are insensitive to the viscosity and lie in the low range of the Raman 
spectrum. Thus, the backbone modulates the motion of base-pairs in accord with the 
environment medium. 

We aim at finding effects of resonance attenuation of hypersound propagating 
in a sample of molecules of the DNA. To that end one may employ the interaction 
between solvent and molecules of the DNA to generate phonon modes of the DNA 
by pumping GHz-excitations in the liquid. The strong viscous interactions of the 
molecules and the solvent may result in a dragging that could promote torsional 
phonon modes generating the interstrand ones, and result in the additional absorption 
of hypersound. Thus, the study of solutions of the DNA could provide important 
information as to the nature of the liquid state, precisely through taking into account 
effects of dissipation in solvent in the GHz range. It is worth noting that hypersound 
acoustics has made considerable progress in recent years and has become a powerful 
method in experimental research, [19], [21]. In the context of the DNA it may turn 
out to be a valuable means for studying hydrodynamic phenomenae. 
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2 Librational dynamics of base-pairs 



In this paper we are concerned with the phonon modes of the DNA. To that end we 
employ the theoretical model worked out in our earlier paper |22| , in which we follow 
the guidelines cast by H.Kappellmann and W.Biem, [23]. For the convenience of the 
reader we recall the main points of paper [22[. In considering the dynamics of the 
DNA one has to take into account: 

1. the DNA having the two strands; 

2. the base-pairs being linked by the hydrogen bonds; 

3. the helical symmetry. 




Figure 1: Base-pairs: adenine — thymine and guanine — cytosine. Dashed lines desig- 
nate hydrogen bonds. 

We shall utilize a one-dimensional lattice model of the DNA which accommodates 
these requirements, in agreement with paper |TB1. We follow the scheme introduced 
by El Hasan and Calladine, [M]. Aiming at a qualitative description of the DNA 
dynamics we use a simplified set of variables, and describe the relative position of 
the bases of a base-pair by means of the vector Y; which is equal to zero when the 
base-pair is at equilibrium, see FIGsdJ [2l The relative positions of different base- 
pairs are described by torsional angles of the sugar-phosphate backbone. They 
correspond to deviations from the standard equilibrium twist of the double helix, so 
that a twist of the DNA molecule, which does not involve inter-strand motion or 
mutual displacements of the bases inside the pairs, is determined by the torsional 
angles (pn of rotation of the base-pairs about the axis of the double-helix. Thus, at 
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equilibrium, a base-pair can be transformed into the next one by rotation through 
the pitch angle and simultaneous translation through a distance of neighbouring 
base-pairs. This picture of the conformation of the DNA is in full agreement with 
that of paper [TH] . 




Figure 2: Angles indicate deviation from equilibrium. Angle f2 = 7r/5 denotes 
the pitch of the double helix. Vectors Nn and Dn designate the unperturbed and 
perturbed states of the bases in a base-pair. Vector y„ equals to the difference between 
Dn and iV„, that is y„ = i3„ - 7V„. 

We suppose that the size of DNA molecule is small enough that it can be visualized 
as a straight double helix, that is not larger than the persistence length. Hence the 
number of base-pairs, N < 150, approximately. Consequently, the twist energy of the 
molecule is given by the equation 



E 



2 Za-^ 



in which /„ is the moment of inertia of the n — th base-pair, and t„ are the twist 
coefBcients, which may change from one base-pair to another. Interstrand motions 
should correspond to the relative motion, or libration in terms of paper |18j . of the 
bases inside the base-pairs, therefore the kinetic energy due to this degree of freedom 
may be cast in the form 

n 

where Mn is the effective mass of the n — th base-pair. 

For each base-pair we have the reference frame in which z-axis corresponds to the 
axis of the double helix, y-axis to the long axis of the base-pair, x-axis perpendicular 
to z- and y- axes (see Fig. 1 of paper [21] )• At equilibrium the change in position of 



4 



neighbouring base-pairs is determined only by the twist angle n of the double helix. 
We shall assume fl = 27r/10 as for the B-form of DNA. To determine the energy due 
to the inter-strand displacements we need to find the strain taking into account the 
helical striicture of our system. Aiming at a simplified picture of inter-strand motion, 
we confine ourselves to the torsional degrees of freedom of the double lattice and 
assume the vectors F„ being parallel to x-y plane, or two-dimensional. Consider the 
displacements determined within the frames of the two consecutive base- 

pairs, n, n-|-l. Since we must compare the two vectors in the same frame, we shall 
rotate the vector l^+i to the frame of the n-th base pair, 

= i?-i(^)y„+i 

Here R^^((p) is the inverse matrix of the rotation of the n-th frame to the (n4-l)-one 
given by the equation 

cos if — sin If 
sin (p cos 



(1) 



The matrix i? is 2 by 2 since the vectors Yn arc effectively two-dimensional. Then the 
strain caused by the displacements of the base-pairs is determined by the difference 

back 



It is important that the angle Lp is given by the twist angle, 0, describing the 
double helix, in conjunction with the torsional angles (^„, so that 

ip = Q. + ipn+l - ipn 

Therefore, the energy due to the inter-strand stress reads 

r ^ ^ 1 2 



E 



2 " ^2a2 



where fc„ are elastic constants determining the librational motion of base-pairs. It 
is important that they may change from one base-pair to another. The equilibrium 
position of the double helix is the twisted one determined by il, with all (pn being 
equal to zero. Combining the formulae given above we may write down the total 
energy of the DNA molecule in the form 



^ = E 



fn+1 — fn) Yn+1 — Y„ 



2 " 



(2) 
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in which fc„ and a arc the torsional elastic constants and the inter-pairs distance, 
correspondingly. In summations given above n is the number of a site corresponding 
to the n-th base-pair, and n = 1,2, . . . ,N, N being the number of pairs in the segment 
of the DNA under consideration. The last term, e„y^/2 accommodates the energy of 
the intcr-strand separation due to the slides of the bases inside the base-pairs. 

Thus, we have the model of a molecule of the DNA which may be visualized 
as (1 + e) — dimensional one, in the sense that it is one dimensional, from the formal 
point of view, and at the same time accomodates librational motions of base-pairs into 
"two transversal dimensions" outside the axis of the duplex. This is due to Yn being 
directed outside the axis of a molecule, while scalar angles describe the twist of the 
sugar-phosphate backbone. Our next step is to split the above (1 + e) — dimensional 
lattice into three interacting linear chains. To that end we shall cast our variables in 
a complex form. Let us introduce complex quantities Zn in accord with the equations 



Figure 3: Lattice of vectors Yn of relative positions of bases in a base-pair. 

— O^n > ) ~^ — Y^ + iY^^ 

Then we may cast the expression for the energy, E = T -\- U, in the form 

JV .2 N . 
rp _ ^ ^nfn _|_ ^ rnnZnZn 

n=0 ^ n=0 ^ 

for the kinetic energy, and 

JV ^ N-l N-l , 

^ = E + E + E ^ 

n— n— n— 

for the potential one. On using the substitution 

— ^ yn 

we may cast the above equations in the form 



- e^(" + ^"+1 - 'fin), 



AT r -9 N 

fi 

2 2 

n=0 n=0 



^^ = E^ + E' 
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and 



N ^ N-1 N-1 , 2 

n— n— n— 



We cosider only small deviations from the equilibrium, and assume that difference 
i^n+i ~ fn) is small. Therefore, we have the equation 



In what follows we shall neglect terms of the second order in (p. 



Within the above approximation the Lagrangian equations of motion have the 
form 

a^InVn = T„-l<y5„_l - (t„_i + r„)(^„ + T„l^„+l (3) 
— kn{anhn+l — an+lhn) + fcn-l(aji-l&ri — an^n-l), 

a^rrindn = - a^^nan + fc„-ia„-i - (fc„_i + fc„)a„ + fc„a„+i (4) 
a^mjjn = - a^e„fe„ + fc„_i6„_i - (fc„_i + fc„)6„ + fc„6„+i (5) 

a„ and 6„ being the real and imaginary parts of ?/„. Thus, we split the equations 
of motion into three parts, and cast the system in a form of three interacting one- 
dimensional lattices of variables a„, 6„, (/?„. This mathematical device is crucial for 
the subsequent analysis of the dynamics. 

It is worth noting that, in real life, a molecule of the DNA is not a perfect duplex, 
so that constants M„, J„,fc„,en are, generally, random sequences. Indeed, owing to 
different sets of base-pairs, the conformational parameters of a molecule of the DNA, 
such as angles determining positions of base-pairs, suffer considerable, upto tens of 
per cent, deviations from the regular positions, [24j . Further, their being of the two 
types, adenine — thymine and guanine — cytosine, results in constants M„,/„,fc„,e„ 
having different values for different pairs. Therefore, equations (jS]) — © describe 
the motion of irregular interacting lattices. One can appreciate this circumstance 
according to ones lights. For one thing, it results in formidable difficulties as far 
as the standard analytical treatment is concerned, for another, it is conducive to a 
phenomenon peculiar to molecules of the DNA in real life, that is the formation of 
localized excitations. In this paper we shall prefer the second option. 

We feel that one can understand the phenomenon of localized excitations within 
the framework of the theory for dynamic excitations in random media worked out. 
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years ago , by I.M.Lifshitz and his colleagues, [2H] — [31] (see also the excellent review 
[33]). According to the Lifshitz theory the dynamics of elastic excitations in random 
media, for example lattices of real crystals or amorphous systems, has two important 
properties. Firstly, there exists forbidden bands in the spectrum of elastic excitations, 
such as phonons. Secondly, there are excitations localized in bounded regions of the 
medium, their characteristic frequencies being located at the edges of the forbidden 
zones. The key observation due to I.M. Lifshitz was that defects of media could serve 
centers for the localization of elastic excitations, see [2S] — [32]. In his first paper 
on the subject he considered defects caused by an atom of isotope that has a mass 
different from the regular one, and showed that the defect could generate localized 
vibrations at its side. Later, the approach was further extended to defects of a more 
general kind and higher dimensions, see [33], for example plains of defects, and also 
boundaries of crystals. 

The study of spectra for disordered, or random, lattices has followed mainly two 
paths. A number of papers employed analytical methods, which generally could 
provide only qualitative description of the spectra, see [33]. It was P.Dean, [35], who 
developed a powerful computer technic for analysis situs of the frequency distribution 
of disordered lattices. The results by P.Dean comply with the Lifshitz theory, and 
provide a powerful insight into a specific structure of spectra. 



3 Hypersound spectroscopy. 

Phonon modes may result in the excessive absorption of hypersound in films and 
solutions of the DNA. In fact, the passage of an acoustic wave may promote the 
transfer of molecules from the equilibrium state to a state in which phonon modes of 
molecules of the DNA are excited. The time delays in this process and its reversal 
should lead to a relaxational dissipation of acoustic energy, and an absorption of hy- 
persound at certain resonance frequencies corresponding to frequencies of the phonon 
modes. Equally important, hypersound irradiation of molecules of the DNA could 
make for generatung phonon modes. Both the experimental data, [3] — [7], [TO] . 
and the theoretical arguments, [T], [T3], [H], indicate that the frequencies of phonons 
of the DNA are in GHz range, and therefore one may expect resonance interaction 
between hypersound waves and phonons of the DNA. The availability of hypersound 
transducers, see [TO], [10], suggests that there are technical means so as to employ 
high resolution acoustical spectroscopy for studying phonon modes of the DNA in 
GHz range. 

The alleged attenuation requires a careful choice of right samples of the DNA for 
studying phonon modes. By now films of the DNA are generally employed to that end. 
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Perhaps, this could make for diminishing the attenuation effects mentioned above. In 
fact, solutions of the DNA are not the best proposition owing to the attenuation 
effects being large in this case. But, any way, it should be very interesting to use 
liquid crystallin phases of the DNA, as was done for inelastic x-ray scattering, see 
paper [33] and references therein. 

The attenuation of phonon modes is closely related to the problem of hydration of 
the DNA. It is alleged that there are two relaxational processes in a hydrated molecule 
of the DNA due to the primary and the second hydration shell with relaxational times 
Ti = 4 X 10~^^ and T2 = 2 x 10~^^ sec. The residence time for a molecule of water 
at grooves of a molecule of the DNA is estimated as 0.2 x 10~^ 0.4 x 10~^ sec, 
[lU] . Phonon dynamics having characteristic times 10^^ 10~^^ sec, there is a 
need for using the generalized hydrodynamics in GHz range, |26j . to estimate, for 
example, the size of dissipation forces acting on a vibrating molecule. Presently, this 
problem is hard to solve, and we have got only the information on the dispersion and 
attenuation of sound waves in this range, required for the theoretical treatment of the 
Mandelstam — Brillouin light scattering. 

But, curious enough, one may expect that the strong attenuation of some phonon 
modes could bring about the generation of others because of the attenuation due 
to the hydration shell mentioned above. The latter primarily involve only external 
regions of the duplex, or according to paper [T5] does not directly affect the librations 
of base-pairs. Hence, we may suggest that an acoustic wave could drag the molecule, 
promote its rotational motion, and generate librations of base-pairs owing to the 
internal dynamics of the duplex. We may infer from equations ([3|) — ([5]) that our 
model allows for the process. The problem is to write down a reasonable force of 
molecular-liquid interaction. 

For the convenience of numerical simulation it is worthwhile to choose appropriate 
scales for length, mass, and time, which agree with the conformational structure of 
the DNA. We shall take the following quantities as the DNA units 

• M = 10~^^ gr as unit of mass, by the order of magnitude close to the mass of 
a base-pair; 

• L = 3 x 10^^ cm as unit of length, close to the distance between neighbouring 
base-pairs; 

• T = 10"^'^ sec as unit of time; corresponding to the upper edge of phonon 
frequencies of the DNA. 
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It is important that the sound velocities, c, of a molecule of the DNA are of 
order 10^ cm/sec, according to various experimental and theoretical estimates, [5], 
[5^ . In the units indicated above we have therefore c ^ 0.3. The moment of inertia, 
/„ , of a base-pair is equal by orders of magnitude m„ , where R is the radius of 
a molecule of the DNA, that is ~ 10 A. Therefore, we have ~ 10 in the DNA 
units introduced above. To obtain numerical estimates for the constants A:„, r„, e„ we 
shall require that the values of sound velocities, c, and frequencies, w, of librations 
inside base-pairs, be of the orders of magnitude c — 0.3 and 1, respectively. Since 
kn / nin ^ and t„ / /„ ^ c^, we have fc„ ^ 0.1 and r ^ 1, respectively. For the 
libration frequency we have e„ / m„ ~ 1, if we assume that it be of order 1 GHz. For 
the viscosity coefficient of water we have, accordingly, fi « 0.3. Summarizing, we have 
the following characteristic quantities expressed in the DNA units 

c ~ 0.3 

~ 10 

- 0.1 

- 1 
~ 1 

fi - 0.3 

As was mentioned above, the hydrodynamical theory, presently, does not provide 
reliable theoretical instruments for studying the interaction between a molecule and 
solvent, and one should turn to a rule of thumb to find a solution. The interaction 
could be small enough, and therefore proportional to the velocity of a molecule. Con- 
fining us to the picture given by paper [18j . we may suggest that the viscous force 
imposed on the n-th cite of the molecule should have the form 

dn = - Jifn (6) 

where 7 is a dissipative constant. The size of 7 is difficult to assess. Estimates based 
on the conventional picture of liquid motion, that is the Navier — Stokes one, are 
apparently far from reality because of small characteristic ("nano") times, and the 
specific structure of water close to the molecule (see [lO]), so that there is a problem 
as to whether it is reasonable to employ the usual hydrodynamic viscosity coefficient 
in this situation. If we assume that it is possible, the dimension analysis gives for the 
moment of viscousity at site n 

or 

^ (X fia^ R^ 



In 
Tn 
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In the DNA units introduced above ii « 0.3, a = 1., i? w 3 so that we have 
7« 10. 

In choosing the drag due to the action of a sound wave on the molecule we have 
to take into account that the wavelength of a hypersound wave in GHz region is by 
several orders of magnitude larger than the molecular segment. The latter being of a 
few hundred A, the force depends only on time. We assume also that the torsion of 
molecule due to the force is small, and may take it in the form 

f ^ AipnCos{ujt) (7) 

The main point is to determine the nature of the force and to assess its size. If we 
confine ourselves to the macro scopical Navier — Stokes hydrodynamics, there are two 
different pictures surfacing. 

First, there could be an effect similar to that of the Rayleigh disc, [36]. To see 
the point let us neglect the dissipation and attribute the force to the streamlines of 
flow round a molecule. There is the turning — moment Q experienced by an oblate 
ellipsoid, or disk, in a vibrating fluid, |57] . 

4 , o 
Q = - pW'v sin 2^9 
3 

The dimension analysis shows that there is the general formula for the turning — 
moment 

Q oc pL^ v'^Lp 

where L is the characteristic length of a body, v velocity and p the density of the 
fluid. In the case we are considering, we may suggest that the rotation angle is small 
so as to employ ip instead of its sinus. Since the coefficient at ip is the product a? /„, 
we obtain the expression 

A oc pR^ (? 

Now we are in a position to assess the value of A as regards the energy pumped in 
a sample. Let us recall that the density Je of the energy current is given by the 
equation Je = pcw^, c being the velocity of sound in liquid. We are interested only 
in rough estimates. Therefore we assume c ~ lO^cm/sec, p ^ \ gramme. For the 
energy current = 1 W jcm? we obtain v — IQcm/sec, or v ^ 0.3, 10^^ in the DNA 
units. In the equation given above R is the radius of the molecule, that is ^ lO^*" cm, 
or i? ^ 3 in the DNA units. Consequently, we have A ^ lO^*". Compared with 
the values of /„, 7, it is too small to produce any appreciable effect. In fact, even 
for Ie ^ I kW/am? we obtain only A ^ lO^'*, which is also too small. Therefore, 
the effect of the Rayleigh disk could not have any bearing on phonon modes of the 
DNA. But it is worth noting that the above estimate is due to the use of the laws 
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of macro scopical hydrodynamics, whereas there are no definite conclusions as to their 
validity on microscopical scale. 

The second option is provided by the viscous interaction of the solvent with the 
molecule. We are looking for an analogue of the familiar Stokes force " dnrjiR^v ". 
Using the dimensional analysis we obtain the expression 

AccfiR^a^v (8) 

For Je = I Wlcrr? it provides the estimate A ~ 10^"', or A ~ 3 • 10^^ for Je — 
1 kW/cm"^. 

The above estimates apparently preclude any opportunity for the observation of 
phonon modes of the DNA by means of hypersound pumping. Similar arguments, 
|13] , are generally put forward as regards phonon modes studied by the sub-millimeter 
absorption spectroscopy. But nonetheless the sub THz-phonon modes are observed, 
see [TT] and the references therein. We feel that this discrepancy between theory 
and experiment is likely to be due to the use of conventional hydrodynamics in the 
region where it does not work properly. In particular, it results in the values of 
the coefficients A and 7 in equations ([3]) — ([5]) that precludes the existence of phonon 
modes. Therefore, in what follows we consider consequencies of A and 7 being outside 
the range prescribed by the Navier — Stokes hydrodynamics, and choose appropriate 
values that permit the existence of phonon modes . 



4 Numerical simulation of phonon modes. 

The equations of motion are obtained by inserting the forces ([6|), ([7]) in equation ([3|). 
The equations are nonlinear and require numerical simulation for their studying. For 
the numerical simulation of the above equations we have used the algorithms of Verlet, 
LeapFrog, and the explicit and implicit Adams algorithms. We have systematically 
made comparisons between results provided by different algorithms so as to escape 
possible mistakes and numerical artifacts. We have taken the integration step 0.01, 
or 0.001 for verification, the unit of time being the period of external force generated 
by sound pumping.. 

We are going to simulate the generation of phonon modes by pumping hypersound. 
According to the equations of motion it is allowed owing to the term given by equation 
(O. Thus, we shall generate firstly a (/3-mode, and the latter shall promote the a„- 
and the 6„-ones, according to equations Q and ([5]). As far as these modes are 
concerned, we have a parametric excitation through the interaction terms in equations 
O,®, and, as is shown below, there may emerge a parametric resonance at a certain 
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Figure 4: Minimal amplitude A against dissipation 7. Resonance frequencies 1.3 
(dashed) and 1.4 (solid). The analytical fit is linear approximations I.2967 + 0,019 
and 1.3957 + 0, 033, respectfully. 

frequency uJex of the excitation pulse. Here it should be noted that the usual theory 
of parametric resonance for the harmonic oscillator, |36| . cannot be directly employed 
because the equations Q — ©of motion for the interacting a„, 6„ chains are 
nonlinear, and we need to find the resonance frequencies by "trial and error" , Lifshitz's 
theory providing at this point general recommendations. 

It is important that there exists the threshold Amin that depends on the dissipative 
coonstant 7, the resonance taking place only for amplitudes of acoustic field which 
are large enough, A > Amin- 

If the frequency of the acoustic wave is very close to the resonance one, Amin 
depends linearly on the dissipative constant 7, see FIGUl ^nd for small 7, Amin is 
also small, so that there is no threshold value at 7 = and Amin = 0. If the excitation 
frequency appreciably deviates from the resonance one, there is a threshold for Amin 
at 7 = 0, see FIGlSl Thus, the numerical simulation of equations ([3]) — ([5]) is in 
agreement with the general theory of parametric resonance. 

Localization of a phonon mode excited by the resonant acoustic wave is illustrated 
in FIGlUfor zero-dissipation case 7 = 0, in FIG|9]for the dissipative constant 7 = 0.01, 
and in FIGlTT]for 7 — 0.02. Our results thus in agreement with Van Zandt and Saxena, 
[17] , who predicted the phonons in the submillimiter range corresponding to localized 
excitations spread over several base pairs. 
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Figure 5: Threshold Amin against the dissipation 7. The excitation frequency equals 
1.55 and deviates from resonance one. Number of sites N — 300, mass m — 
1; elastic constants r ~ 1 or 1.5 (equal probabilities); e = 1; k = 
0.1 or 0.15 (equal probabilities); / = 10. The value of Amin at 7 = is not 
equal to zero. 



Figure 6: Ranges of the resonance frequencies corresponding to the excitation ampli- 
tudes 0.01, 0.02, 0.03, 0.04, 0.05 (thick horizontal segments); e = 1; fc = 0.1 or 0.15 
(equal probabilities); / = 10; 7 = 0. 
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5 Conclusion: a need for acoustic spectroscopy of 
the DNA 

We feel that the acoustic spectroscopy could be helpful in analyzing the conforma- 
tional structure of the DNA through the analysis of its phonon modes. The model 
we have constructed qualitatively agree with the experimental data by providing the 
existence of the phonon modes and the reasonable orders of magnitude for their fre- 
quencies. It implies a phonon localization due to the random structure of the duplex, 
and shows that it would be worthwhile to study the action of hypersound on molecules 
of the DNA. It is important that absorption peaks for hypersound could be expected 
at frequencies corresponding to the parametric resonance of phonon modes of the 
DNA under hypersound pumping. In this respect the parametric resonance could be 
a valuable means for studying the vibrational modes of the DNA. 

There still remains a problem of accommodating the dissipative effects. The stan- 
dard Navier — Stokes model promotes a size of dissipation that exclude the existence 
of phonon modes, but they are observed in real life, [TT]. The circumstance challenges 
theoreticians to working out an adeqate framework for the hydrodynamics at molec- 
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Figure 8: Dissipation 7 = 0. The DNA units are employed. Localized excitation on 
a random lattice; N = 300, to = 1;t = 1 or 1.5 (equal probabilities); e = l;k = 
0.1 or 0.15 (equal probabilities); I = 10; amplitude A = 0.08; excitation frequency 
Lj = 1.5. Time elapsed 300 periods of the acoustic wave. The graphs correspond to 
the a—,b—, modes (view from the top). 
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Figure 9: Dissipation 7 = 0.01. The DNA units are employed. Localized excitation 
on a random lattice; N — 300, m = 1: r = 1 or 1.5 (equal probabilities); e = l;k = 
0.1 or 0.15 (equal probabilities); / = 10; amplitude A = 0.1; excitation frequency 
IV = 1.489. Time elapsed 1000 periods of the acoustic wave. The graphs correspond 
to the a—, b—, (^— modes (view from the top). 
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Figure 10: Dissipation 7 = 0.02. The DNA units are employed. Localized excitation 
on a random lattice; N = 300, m = 1; r = 1 or 1.5 (equal probabilities); e = 1; fc = 
0.1 or 0.15 (equal probabilities); / = 10; amplitude A = 0.1; excitation frequency 
ui = 1.489. Time elapsed 1300 periods of the acoustic wave. The graphs correspond 
to the o— , b—, modes (view from the top). 
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Figure 11: Dissipation 7 = 1. The DNA units are employed. Localized excitation on 
a random lattice; TV = 300, m = 1;t = 1 or 1.5 (equal probabilities); e = l;fc = 
0.1 or 0.15 (equal probabilities); / = 10; amplitude A = 1.7; excitation frequency 
LJ = 1.489. Time elapsed 400 periods of the acoustic wave. The graphs correspond to 
the a—,b—, modes (view from the top). 
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Figure 12: Dissipation 7 = 0.8. The DNA units are employed. Localized excitation 
on a random lattice; N = 300, m = 1; r = 1 or 1.5 (equal probabilities); e = 1; fc = 
0.1 or 0.15 (equal probabilities); / 10; amplitude A — 1.3; excitation frequency 
oj = 1.489. Time elapsed 760 periods of the acoustic wave. The graphs correspond to 
the a—,b—, modes (view from the top). 
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ular scale. The part played by the dissipation, perhaps different from the Navier — 
Stokes one, could be crucial, and torsional modes of the backbone could make for 
understanding its nature on the nano-scale. Hypersound acoustic spectroscopy could 
be instrumental in this respect. To that end liquid samples may turn out to be more 
interesting, solutions of the DNA providing a probe into the physics of liquid state. 
The observation of the parametric resonance we have indicated could serve a means 
for understanding the hydrodynamics on molecular scale. 

This paper was partially supported by the Russian Foundation for Fundamental 
Research, grant # 09-02-00551-a. 
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